N PS  ARCHIVE 

1966 
RADECKI,  R. 


SHALLC  .OACOUSTTC  CO 


^ 


RICHARD  A  ICKI 


MAS:  ,TTS  INSTITUTE  OF  TECHNOLOGY 


MAY,   1966 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  93943*101 


iate  School 


SHALLOW   WATER  HYDROACGUSTIC   COMMUNICATIONS 

by 


Richard  Anthony  Radecki 

// 

B.S.,   United  States   Naval  Academy 

1959 


SUBMITTED   IN  PARTIAL  FULFILLMENT   OF  THE  REQUIREMENTS 
FOR  THE  DEGREES    OF  NAVAL  ENGINEER  AND   MASTER 
OF  SCIENCE.    (BIBGTBICAL  ENGINEERING) 

at   the 

MASSACHUSETTS    INSTITUTE   OF  TECHNOLOGY 
June ,    1966 


Signature   of  Author    

Department   of  Naval  Architecture  and   Marine  Engineering, 

May  20,    I966 


Certified  by   

Thesis  Supervisor 


Accepted  by    

Chairman,   Departmental  Committee   on  Graduate  Students 


,duate  School 
rata 


DUDLEY  KNOX  LIBRARY 

NAVAL  POSTGRADUATE  SCHOOL 

MONTEREY,  CA  93943-5101 


SHALLOW  WATER  HYDROACCUSTIC  COMMUNICATIONS 

by 
RICHARD  ANTHONY  RADECKI 


Submitted  to  the  Department  of  Naval  Architecture  and 
Marine  Engineering  on  May  20,  1966  in  partial  fulfillment 
of  the  requirements  for  the  Master  of  Science  degree  in 
Electrical  Engineering  and  the  Professional  degree, 
Naval  Engineer. 

ABSTRACT 

Hydroacoustic  communication  has  not  been  used  ex- 
tensively at  ranges  in  excess  of  ten  miles.   A  study  of 
tie  shallow-water  and  surface  ducting  channels  indicate 
that  reliable  sixty  word  per  minute  communication  at  ranges 
of  fifty  miles  is  possible  under  a  variety  of  conditions, 
with  reasonable  power  levels  and  equipment  complexity.   The 
channel  can  be  approximated  by  a  Rayleigh  fading  signal  with 
additive  Gaussian  noise.   Five  frequency  multiplexed  binary 
channels  are  compared  with  a  fifteen  channel  system  utilizing 
frequency  diversity,  and  with  a  fifteen  channel  system  util- 
izing coding.  A  method  of  calculating  the  diversity  error 
for  components  with  unequal  signal  and  or  noise  levels  is 
worked  out.   Frequencies  must  be  above  900  cps  if  shipping 
noise  is  significant,  and  below  1500  cps  to  minimize  at- 
tenuation,  ^o  optimize  several  transmitter-receiver  para- 
meters, such  as  the  signal  interval,  more  long  range  ex- 
perimental data  is  required. 
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INTRODUCTION 

Except  for  a  few  very  short  range  or  very  Iovj  rate 
systems,  the  communications  capability  of  ocean  acoustics  has 
not  been  exploited.   Voice  modulation  has  been  used  extensively 
for  ranges  of  a  few  miles  and  CW  at  slightly  longer  ranges. 
Single  explosive  charge  signals  have  been  sent  over  a  thousand 
miles  via  deep  ocean  channels.   Aside  from  sonar  these  few 
examples  cover  most  of  the  communications  development.   With 
the  increasing  interest  in  nuclear  submarines,  oceanography, 
and  the  industrialization  of  underwater  real  estate,  the  use- 
fulness of  a  reliable  system  which  does  not  depend  on  radio 
or  cables  is  apparent. 

In  this  thesis  a  single  such  problem  will  be  con- 
sidered, that  is,  communication  between  two  stations  in  rela- 
tively shallow  water  using  transducers  near  the  surface.   A 
range  of  fifty  miles  and  a  data  rate  of  sixty  words  per  minute 
will  be  used  to  evaluate  various  systems.   If  such  a  system 
were  feasible  it  x^ould  provide  a  good  link  between  ships  ope- 
rating in  the  same  geographic  area,  or  between  an  oceano^raphic 
vessel  and  its  remote  data  collecting  stations,  or  any  number 
of  other  purposes. 

Part  I  of  this  thesis  will  be  a  simplified  description 
of  the  shallow  water  acoustic  channel,  which  will  point  out 
the  parameters  important  to  a  communications  system.   These 
include  ambient  noise  levels,  transmission  losses,  and  channel 
stability.   xhese  parameters  will  be  used  to  arrive  at  a 


random  process  model  of  the  channel. 

In  Part  II  this  model  will  be  used  to  develop  a  com- 
munication system.   First  a  single  binary  channel  will  be 
developed.   This  will  then  be  extended  as  required  to  obtain 
the  desired  data  rate  at  a  range  of  fifty  miles  with  an  ac- 
ceptable error  rate,   ^ince  there  are  many  ways  to  handle  the 
problem  only  one  or  two  representing  current  technology  will 
be  discussed. 

In  Part  III  the  feasibility  of  hydroacoustic  communica- 
tions under  various  ocean  conditions  will  be  evaluated.   The 
assumptions  used  in  the  first  two  parts  and  their  effects  will 
also  be  evaluated  at  this  point. 

Part  IV  will  summarize  the  results  and  indicate  areas 
requiring  additional  research  and  or  development. 


Part   I 
SHALLOW  -WATER  EYDROACOUSTICS 

1.  General.   In  shallow  water  there  are  two  principle  methods 

of  acoustic  transmission,   The  first  type  is  the  surface  channel, 
In  this  channel  a  large  portion  of  the  acoustic  energy  is 
trapped  in  a  duct  near  the  surface.   Transmission  in  this 
duct  is  by  surface  reflection  and  upward  refraction.   The  second 
type  is  transmission  by  both  surface  and  bottom  reflections. 
The  particular  type  of  transmission  which  villi  occur  is  pri- 
marily dependent  on  the  vertical  profile  of  the  velocity  of 
sound  in  the  ocean  area.  The  velocity  of  sound  in  sea  water 
varies  primarily  with  temperature  and  to  a  lesser  decree  with 
pressure  and  salinity.   Over  long  ranges  it  is  reasonable  that 
both  types  of  transmission  may  occur  for  part  of  the  distance. 

2.  The  Surface  Channel,   ^he  density  of  water  decreases  with 
temperature.   r^his  means  that  in  general  the  water  near  the 
surface  is  warmer.   However,  if  a  cooling  condition  exists, 
the  surface  water  may  be  cooler  or  at  the  same  temperature 

as  deeper  water.   This  inversion  causes  the  surface  channel  to 
exist.   The  velocity  of  sound  increases  with  increasing  pres- 
sure and  with  increasing  temperature  (1).   Therefore,  an 
inversion  or  isothermal  layer  will  have  a  positive  velocity 
gradient,  i.e.,  velocity  of  sound  increases  with  depth.   Ac- 
cording to  Snell's  law,  acoustic  waves  in  this  layer  will  be 
bent  upward,  reflected  from  the  surface  and  bent  upward  again. 
Figure  1  is  a  diagram  of  this  type  of  channel.   It  has  been 
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Figure  1.   Surface  ducting  channel. 
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Figure  2.   Spherical  and  cylindrical  spreading  losses. 


shown  and  experimentally  verified,  that  a  minimum  layer  depth 
is  required  to  trap  a  given  frequency  signal  (2).  This  depth 
is: 

Lmin  =  2.?  x  103  (fg)_1/3 

f  =  frequency  (cps ) 

g  =  duct  velocity  gradient  (sec.-1) 
For  a  1000  cps  signal,  and  a  0.020  meter/second/meter  velocity 
gradient  this  minimum  layer  depth  is  30/f  meters. 

There  are  several  types  of  transmission  losses  which 
must  be  considered  in  this  channel.   These  are:   spreading 
losses,  reflection  losses,  and  absorption  losses.   The  spread- 
ing loss  will  be  broken  into  spherical  and  cylindrical  com- 
ponents.  The  energy  leaves  the  source  initially  as  a  spherical 
wave.   Beyond  a  certain  point  the  energy  is  essentially  trapped 
in  the  layer  and  the  spreading  loss  becomes  cylindrical.   For 
purposes  of  calculation  the  transition  point  will  be  taken  as 
half  the  horizontal  range  of  a  ray  leaving  the  source  at  the 
layer  depth  horizontally  and  just  reaching  the  surface,   '■'•his 
is  shown  in  figure  2. 

The  second  type  of  loss  to  be  considered  is  the  surface 
reflection  loss.   The  reflection  coefficient  for  the  water- 
air  interface  is  approximately  minus  one.   This  is  true  for 
calm  water,  but  waves  cause  a  certain  amount  of  scatter  loss. 
This  loss  coefficient  has  been  determined  experimentally  (3)» 
The  following  formula  describes  the  loss  function  for  a  fre- 
quency-waveheight  product  less  than  J: 


<*,=  -10  log/o  (1-0.  1^9  (fH)3/2)   db  per  bounce 

f  =  frequency  in  kc/sec. 
H  =  waveheight  in  meters 

The  third  type  of  loss  to  be  considered  is  absorption. 
This  takes  into  account  primarily  the  thermal  and  viscous 
forces  which  are  otherwise  neglected  in  deriving  the  wave 
equation  of  acoustics  in  a  liquid  medium.   For  the  ranges  of 
interest  this  loss  is  negligible  belovi  frequencies  of  2  kc. 
A  suitable  formula  for  absorption  losses  which  is  accurate 
from  frequencies  of  2-25  kc  is  ('+)  : 

°<e=  720Q  *  10~  f  +  3.35  x  10~8  f2""nepers  per  meter 
7750  +  f2 

this  gives  a  transmission  loss  of: 

L  =7.2  80  x  103x  <*2   db/km. 

If  it  is  assumed  that  the  velocity  of  sound  does  not 
vary  in  the  channel  except  linearly  with  depth,  some  rather 
simple  equations  can  be  derived  which  specify  the  geometry 
of  the  path  in  depth  and  range.   Hef erring  to  figure  3. 
these  formulas  are  (5): 

d  =  transducer  depth  (meters ) 

L  =  layer  depth  (meters) 

R  =  horizontal  range  (meters ) 

a  =  velocity  gradient  (sec."1) 

^__  900  _  angle  of  incidence 

Y  =  vertical  depth  used  for  calculation  (meters) 

c0=  lowest  velocity  of  sound  in  the  layer  (meters/sec.) 
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Figure  3.   Surface  channel  ray  path  calculation  geometry. 
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Figure  4.   The  SOFAR  channel. 
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If  we  now  let  Y  =  L  (transducer  at  the  surface),  and  -e-.  =  0 

then  2R  is  th ©maximum  distance  between  reflections,  and  the 

minimum  number  of  reflections  over  a  distance  D  is: 

minimum  number  of  reflections  =  D_ 

2R 

Since  the  water  depth  does  not  enter  any  calculations 

for  the  surface  channel,  it  is  also  a  suitable  deep  water  model 

as  long  as  it  is  uninterrupted. 

3.   The  Bottom  Reflection  Channel.   The  bottom  reflection 
channel  only  exists  in  shallow  water.   In  very  deep  water  the 
temperature  reaches  a  minimum  at  some  intermediate  depth  and  is 
isothermal  down  to  the  bottom.   This  gives  rise  to  a  minimum 
velocity  of  sound,  which  generally  occurs  at  depths  of  about 
two  thousand  feet.  Sound  waves  are  reflected  upward  below 
this  depth,  and  downward  above  this  depth.   %is  is  the  SCFA.R 
channel  shown  in  figure  4.   It  will  not  be  considered  in  this 
thesis.   For  the  shallow  water  channel,  ocean  depths  up  to 
600  meters  will  be  considered.   This  includes  about  ten  per 
cent  of  the  oceans,  mostly  along  the  continental  shelf  (6). 

In  shallow  water  the  SOFAR  channel  rarely  exists,  and 
the  effects  of  both  bottom  and  surface  reflections  are  signi- 
ficant.  Some  typical  linearized  velocity  profiles  are  shown 


X 


in  figure  5.    For  uniform  velocity  gradients  the  ray  path 
equations  as  described  for^surface  channel  may  be  used  for 
the  bottom  reflection  channel  path  calculations.   For  a,  b, 
and  c  of  figure  5,  these  calculations  would  have  to  be  made 
in  several  steps. 

Except  for  the  addition  of  a  bottom  reflection  loss 
the  bottom  reflection  channel  losses  are  the  same  as  those  of 
the  surface  channel.   Because  of  the  greater  channel  depth  it 
is  obvious  that  the  spherical  spreading  loss  will  extend  for  a 

tor 

greater  distance,   ^-his  will  be  computed  as  R  vice  R/2  as^the 
surface  channel.   In  general,  experimental  data  shows  that 
the  total  losses  for  this  type  of  channel  approximate  spherical 
losses  for  the  whole  transmission  distance.   Aie  surface  channel 
has  total  losses  somewhat  lower  than  spherical  spreading. 
The  bottom  reflection  loss  is  a  complex  function  of  the  bottom 
material  and  its  structure.   For  grazing  angles  less  than 
fifteen  degrees  however,  the  losses  are  small.   A  good  expe- 
rimental figure  is  0.3  db  per  bounce  (7).  For  larger  grazing 
angles  the  loss  becomes  quite  large.   Figure  6  is  a  transmis- 
sion path  diagram  of  the  bottom  reflection  channel.   Note  that 
rays  which  graze  the  bottom  at  about  15  degrees,  leave  the 
source  at  a  much  smaller  angle. 

Sample  calculations  are  made  for  both  the  surface 
channel  and  the  bottom  reflection  channel  in  Appendix  A. 
Several  computations  are  plotted  in  figure  7  together  with 
some  experimental  data  for  comparison.   In  all  cases  the 
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Figure  5.   Shallow  water  velocity  profiles. 
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Figure  6.  The  bottom  reflection  channel 
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Figure  7.   Sound  pressure  level  at  50  miles  as  a  function 
of  frequency.   All  data  corrected  to  a  source  level 
•  of  179  db.  at  a  radius  of  one  meter.   References  for 
experimental  data  in  parentheses  . 
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sources  were  corrected  to  a  sound  pressure  level  of  185  db 
(ref.  0.00002  newt./m2)  at  a  lA  meter  radius. 


SPL  =  20  log1Q  P 

Pc 


o 

This  is  low  enough  to  permit  simultaneous  transmission  over 
each  of  ten  channels. 

k.      Noise .   Ambient  sea  noise  is  a  complex  array  of  processes 
which  are  not  all  fully  understood  (13).   Most  noise  sources 
are  very  frequency  dependent.   Below  1000  cps  ship  and  traffic 
noise  is  the  most  significant.   It  usually  peaks  between  300 
and  800  cps  and  falls  off  rapidly  at  higher  frequencies.   For 
all  frequencies  of  interest  above  1000  cps  the  ambient  noise 
is  primarily  dependent  on  sea  state  for  all  but  calm  weather. 
The  primary  sources  of  this  noise  are  thought  to  be  wind,  spray, 
cavitation,  and  precipitation.   Other  sources  are  marine  life, 
molecular  and  thermal  agitation  and  earthquakes.   Table  1 
gives  the  average  relations  between  wind  force,  wind  velocity, 
wave  height,  and  sea  state,   figure  8  is  a  plot  of  the  average 
ambient  noise  in  the  ocean  versus  sea  state.   The  highest  ex- 
pected values  are  about  3  db  above  the  average.   A  possible 
explanation  is  an  irregular  gusting  wind.   It  should  also  be 
noted  that  deep  water  noise  levels  are  generally  5  db  lower 
than  shallow  water  levels. 

If  a  transducer  is  to  be  carried  on  a  ship,  the  ship's 
self  noise  is  very  important  (1^).   At  low  speeds  this  will 
correspond  to  the  ambient  noise  at  sea  state  2,  but  it  may 
contain  some  high  level  discrete  frequencies.   At  speeds  above 
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Figure  8.  Ambient  noise  in  the  ocean(13).  Specific 


shallow  water  data  indicated  as  follows: 
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about  fifteen  knots,  cavitation  noise  can  become  much  greater 
than  the  ambient  noise. 
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Table  1 

Beaufort  Wind  Speed    Wave  Height 
wind      knots  feet 

force    (meters /sec . )    (meters ) 
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5.   Source  Power.   The  output  power  of  a  transducer  is  limited 
by  cavitation  and  acoustic  non-linearity.  Cavitation  occurs 
when  the  local  pressure  at  or  near  the  transducer  is  less 
than  the  vapor  pressure  of  the  water.   A  sample  calculation 
in  Appendix  B  shows  that  this  limits  the  transducer  to  a 
sound  pressure  level  of  195  <ib  (ref.  0.00002  newt./m2)  at  a 
transducer  depth  of  twenty  feet. 

The  source  level  of  the  transducer  is  limited  by  the 
approximations  required  to  obtain  the  acoustic  wave  equations 
from  the  equations  of  fluid  dynamics,  '^he  governing  ap- 
proximation comes  from  the  following  total  derivative  (17 ): 

du  _Ju    .   u^u 
dt  it  ir 


>y 


15 

u  =  fluid  velocity  vector 

t  =  time 

r  =  fluid  particle  position  vector 
If  the  second  term  on  the  right  can  be  neglected,  then  the 
wave  equation  can  be  satisfied.   This  is  worked  out  in  Ap- 
pendix B  for  spherical  waves.   The  result  is  that  the  sound 
pressure  level  at  one  meter  radius  must  be  less  than  2$h   db 
(ref.  0.00002  newt./rr,2). 

6.   Random  Process  Channel  hodel.   An  input-output  relation 
for  the  hydroacoustic  channel  will  now  be  developed.   The 
parameters  which  will  be  considered  are:  signal  energy,  noise 
energy,  frequency  dispersion,  time  dispersion,  and  multipath 
reception. 

Since  the  ambient  sea  noise  is  due  to  a  large  number 
of  small  sources,  it  will  be  assumed  to  be  a  zero  mean, 
Gaussian  random  process  with  spectral  density  Nf .   Nf  is  then 
the  noise  energy  density.   This  noise  is  independent  of  the 
message  sent.  Therefore  the  received  signal  will  be  a 
message  filtered  by  the  channel  plus  additive  Gaussian  noise. 

This  channel  filter  function  is  unknown  however.   To 
arrive  at  a  channel  description,  two  assumptions  will  be 
made.   First,  along  any  given  transmission  path,  the  channel 
is  a  pure  attenuator.  Second,  the  receiver  receives  the 
signal  over  a  large  number  of  paths,  with  the  received  phase 
for  each  path  being  a  uniformly  distributed  random  variable. 
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Wozencraft  and  Jacobs  show  that  this  description  leads  to 
a  random  process  channel  model,  with  the  following  input- 
output   relationships    (18): 

s(t)   =  m(t)  VF' cos  cjQt  input 

r(t)   =  a(t)m(t)  V27 cos    (w0t-£)+n(t )  putput 

in  which 

p(a)   =  2ae-a2  a>  0 

p(4)   =   1  0  ^3^2-rr 

27T 

The  amplitude  a(t)  is  thus  a  Rayleigh-distributed  random 
variable.   Experimental  data  verifies  this  model  in  many 
instances,  but  other  distributions  have  also  been  observed 
(19,20).   For  the  work  to  follow,  a(t)  will  be  assumed  to  be 
Rayleigh-distributed. 

Before  a  communications  system  can  be  developed, 
frequency  and  time  dispersion  must  be  discussed,   frequency 
dispersion  will  limit  the  banid  length..  Experimental  data 
shows  that  in  general,  for  ranges  of  interest,  frequency 
dispersion  is  less  than  2  cps  (21).   It  is  expected  that  the 
system  will  use  a  baud  corresponding  to  a  much  larger  band- 
width than  this. 

Time  dispersion  of  a  signal  is  the  range  of  delays 
observed  over  paths  of  different  length.   It  limits  the  rate 
at  which  signals  may  be  sent  over  conventional  systems, 
■'■'here  is  very  little  experimental  data  on  shallow  water  time 
dispersion.   The  calculations  in  Appendix  C  show  that  about 
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1/5  second  is  the  longest  time  that  need  to  be  considered, 
at  a  range  of  fifty  miles.   This  is  about  ^0  percent  of  the 
dispersion  for  the  SOFAR  channel  which  is  severely  limited 
in  this  respect  (22),  and  it  corresponds  to  an  experimental 
value  observed  in  the  Florida  Straits  (23). 
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Part  II 
C  OMMU N IC  AT I ONS   S  YS  TE M 

1.   Frequency  Selection.   The  optimum  frequency  is  the  one 
which  alloitfs  the  greatest  bandwidth  and  the  signal  to  noise 
ratio.   Neglecting  the  noise,  low  frequencies  are  the  best. 
Spreading  losses  are  about  the  same  for  all  frequencies. 
■For  a  state  3  sea,  surface  reflection  losses  are  0.5  db./ 
bounce  at  800  cps,  0.95  db. /bounce  at  1000  cps,  and  2.0  db./ 
bounce  at  1500  cps.   This  becomes  worse  for  higher  frequencies 
or  higher  seas.   Attenuation  also  increases  with  frequency 
and  is  first  significant  for  the  fifty  mile  range  at  about 
2000  cps.  Therefore,  primarily  because  of  surface  losses, 
an  upper  limit  of  1.5  kc  will  be  used,  A  lower  limit  of  800 
cps  will  be  used.   At  low  frequencies  it  is  possible  to  get 
strong  harmonics  from  shipboard  machinery,   '■''his  is  general- 
ly worst  at  about  500  cps.   One  prime  example  of  low  frequency 
noise  is  propeller  harmonics  which  may  be  very  strong. 

A  plot  of  noise  and  fifty  mile  signal  strength  is 
shown  in  figure  9  as  a  function  of  frequency,   '^his  plot 
uses  the  same  data  as  figures  7  and  8.  A  noise  level  cor- 
rected to  a  10  cps  bandwidth  is  also  shown.   £his  is  suf- 
ficient for  the  correlation  time  of  1/6  second  which  will 
tentatively  be  used.   Experiments  show  that  at/long  ranges, 
the  signal  will  remain  about  80$  coherent  by  autocorrelation 
for  up  to  lA  second  (23). 
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Figure  9.   Expected  signal-noise  separation  in  a  state 
3  sea.   (See  figures  7  and  8.) 
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2»   Binary  Signalling.   A  simple  binary  communications 
system  will  now  be  developed  for  the  shallow  water  channel. 
The  signal  set  will  be  an  orthogonal  pair  of  cosine  waves, 
as  f ollows : 

s0Ct)  =  V?  cos  CJ0t 

si(t )  =  V21  cos  U   t 
This  set  is  considered  the  easiest  to  use  for  this  system. 
Signals  which  are  time  orthogonal  are  not  well  suited  to  this 
channel  since  the  time  of  arrival  is  unknown.   Since  both 
phase  and  amplitude  are  random  processes,  neither  phase  nor 
amplitude  modulation  was  considered  suitable. 

when  the  signal  s0(t)  is  sent,  the  received  signal  is 

of  the  form: 

r(t)  =  Z_  ai  cos^ot+Qi)  +  n(t) 
i=l 

where  the  i*s  represent  the  different  multipaths.   As  noted 

in  Part  I,  this  leads  to  the  following  random  process  model: 

r(t)  =  a(t)  cos(W0t  +  0)  +  n(t) 

in  which  a  and  0  are  random  variables  with  the  following 

probability  density  functions  at  any  instant  of  time: 

2 
p(a)  =  2a  e_a  0  ^  a 

p(0)  =  1/27/  0^0^27/ 

This  model  and  the  following  assumptions  will  be  used  to 
arrive  at  a  suitable  binary  receiver: 

a.  a(t)  and  0(t)  are  constant  over  the  signal  inter- 
val T,  but  vary  randomly  from  interval  to  interval. 

b.  n(t)  is  a  white  Gaussian  noise  random  process, 
of  spectral  density  No,  and  it  is  independent  of  the  signal. 
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c.  The  two  signals,  s0(t)  and  s1(t)  have  equal 
energy  and  are  equi probable. 

d.  There  is  no  intersymbol  interference. 

A  maximum  likelihood  receiver  will  be  used.   For  this 
receiver,  that  message  is  selected  for  which  the  probability 
of  the  received  signal  conditioned  by  the  messages  is  highest. 
The  details  of  this  receiver  have  been  worked  out  by  Turin 
(26),  and  for  more  general  time-varying  cases  "oy   Kailath  (27) 
and  Price  (28).   In  the  case  of  the  narrow  band  signals  chosen 
for  the  signal  set,  the  receiver  was  determined  to  be  a  de- 
vice which  correlates,  squares,  and  sums  the  quadrature  com- 
ponents before  making  a  decision.   The  correlation  receiver 
is  shown  schematically  in  figure  10.   Figure  11  is  an  envelope 
detector  receiver  which  can  be  shown  to  be  essentially  the 
same. 

The  correlation  receiver  can  be  explained  heuristic- 
ally  by  a  geometric  argument.   The  quadrature  noise  com- 
ponents of  the  received  waveform  are  zero  mean  Gaussian  ran- 
dom variables  of  variance  N0/2,  where  N0  is  the  mean  square 
noise  power.  The  quadrature  signal  components  of  the  received 
waveform  are  also  zero  mean  Gaussian  random  variables,  but 
with  a  variance  of  E0/2,  where  E0  is  the  mean  square  signal 
power.   If  a  signal  is  present  the  received  quadrature  com- 
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Figure   10.      Binary   correlation  receiver. 
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Figure  11.   Binary  envelope  detector  receiver. 
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ponents  are  therefore  zero  mean  Gaussian  random  variables 

with  variance  S,*  N.  .   These  Gaussian  random  variables  are 

2 
illustrated  in  figure  12,  in  the  s„ -plane  and  the  s,  -plane 

for  the  case  when  s0(t)  is  sent.   In  each  case  the  amplitude 

of  the  resultant  vector  is  a  Ray lei gh  distributed  random 

variable . 

The  error  probabilities  for  the  likelihood  receiver  are 

developed  in  Wozencraft  and  Jacobs  (29),  by  computing  it 

first  for  the  random  phase  case  only,  to  obtain: 

-E/2N 

P  (£)  =  I  e 

Averaging  this  error  over  a  Rayleigh  amplitude  distribution, 
the  following  error  probability  is  obtained  for  the  random 
phase  and  random  amplitude  channel: 

p(^)  ~  2+  E/N 

Using  this  formula,  the  efficiency  of  the  simple  binary 

channel  may  now  be  determined.  Referring  to  figure  9,  with 

a  10  cps  bandwidth  and  the  weak  signal  line,  table  2  summarizes 

the  binary  channel  error  probability  at  various  frequencies. 

Table  2 

Frequency      E/N  (db. )        E/N  P(£) 

800  29  790  0.0013 

1000  23  200  0.0050 

1250  15  31.5  0.030 

1500  8            6.3  0.12 

3.   Channel  Capacity  and  Coding.   In  1948,  C.E.  Shannon 

showed  that  whenever  signals  are  transmitted  at  some  rate, 

R,  less  than  channel  capacity,  C,  then  an  arbitrarily  small 

error  rate  can  be  achieved  by  coding  with  a  sufficiently 
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Figure  12.   Received  signal  in  the  s0 -plane  and  s,  -plane 
when  message  s0(t)  is  sent. 
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Figure  12A.   The  binary  symmetric  channel 
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Ion/  block  length  (30).   For  the  binary  symmetric  channel 

with  crossover  probability  p,  as  shown  in  f igure  12A,  the 

channel  capacity  has  the  rather  simple  form  (31): 

C  =  l/T   1  +  plog2p  +  (1-p)  log2(l-p)J   bits  per  second 

where  T  is  the  signal  time  required  to  achieve  p.   It  can 

also  be  shown  that  for  n  arbitrary  channels  being  used  in 

parallel,  the  capacity  is: 

n 

c  =   £Ci 
i=i 

the 

whereACi  are  the  individual  channel  capacities.   This  last 
formula  holds  for  simple  parallel  use  in  which  a  decision 
is  made  on  each  channel.   If  the  decision  is  made  on  the  col- 
lective output  probabilities  of  the  individual  channels,  then 
the  channel  capacity  is  somewhat  greater,  since  in  the  first 
method  some  information  is  lost  by  the  intermediate  decision. 

A  sixty  word  per  minute  data  rate  requires  approxima- 
tely thirty  bits  of  information  per  second.   If  an  average 
channel  capacity  of  five  bits  per  second  can  be  obtained,  then 
a  minimum  of  six  channels  are  required  for  effective  trans- 
mission.  At  transmission  rates  very  close  to  capacity  how- 
ever, very  elaborate  coding  schemes  are  required  to  achieve 
a  reliable  transmission.   In  general,  it  can  be  shown  that 
for  proper  coding  the  probability  of  error  is  bounded  by 
functions  of  the  transmission  rate  (32): 

e-n(EL(R)  +  0(n)    ^   p(£)ie"nE(B) 
where  Et(R)  and  E(R)  depend  on  the  channel,  and  n  is  the  code 
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block  length.  For  rates  less  than  channel  capacity,  both 
Ej^R)  and  E(R)  are  positive.  These  functions  will  not  be 
evaluated  for  this  channel. 

For  this  fading  channel,  the  rate  is  limited  by  the 
number  of  parallel  channels.   The  signal  interval  cannot  be 
shortened  to  much  less  than  one  sixth  of  a  second  because  of 
intersymbol  interference.   The  number  of  channels  is  limited 
by  the  useable  frequency  band.   Above  about  15 00  cps,  very 
little  capacity  is  gained  for  the  addition  of  an  extra  channel. 
Assuming  a  signal  interval  of  one  sixth  of  a  second,  and  trans- 
mitting five  information  bits  during  each  interval,  the  ef- 
fect on  the  rate  as  a  function  of  capacity  of  adding  more 
channels  is  summarized  in  table  3t  and  figure  13.   A  minimum 
useable  frequency  of  900  cps  and  a  frequency  separation  of 
20  cps  was  assumed. 

Table  3* 


Channel  # 

frequencies 

?(0, 

ac 

C 

R/C 

1 

900  &   920 

.0028 

5.8 

2 

940  &   960 

.0040 

5.8 

3 

980  &  1000 

.0057 

5.7 

4 

1020  &  1040 

.0076 

5.6 

5 

1060  &1080 

.0100 

5-5 

28.4 

1.06 

6 

1100  &  1120 

.0124 

5.4 

33.8 

.89 

7 

1140  &  1160 

.0154 

5.3 

39.1 

.77 

8 

1180  &  1200 

.0213 

5.1 

44.2 

.68 

9 

1220  &  1240 

.0263 

4.9 

49.1 

.61 

10 

1260  &   1280 

.0333 

4.7 

53.8 

.56 

11 

1300  &   1320 

.0408 

4.5 

58.3 

.51 

12 

1340  &  1360 

.0562 

4.1 

62.4 

.48 

13 

1380  &  1400 

.0685 

3.8 

66.2 

.45 

14 

1420  &  1440 

.0918 

3.3 

69.5 

.43 

15 

1460  &  1480 

.114 

2.9 

72.4 

.41 

*  Based  on  a  power  level  allowing  ten  channels.   For  less 
channels  P(£)  is  better,  for  more  worse  (See  Part  I). 
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Figure  13.  The  effect  of  adding  additional  higher  frequency 
channels  on  capacity  and  rate.   (Transmission  at  30 
bits  per  second.) 


28 


There  are  several  interesting  methods  of  coding  and 
decoding  which  will  be  discussed  after  developing  a  com- 
munications system  for  the  hydroacoustic  channel.   In  the 
next  section,  the  system  will  be  developed  first  using  simple 
transmission  and  then  diversity  and  block  coding.   Block 
coding  is  taking  a  block  of  input  digits  and  coding  them  into 
a  longer  group  of  coded  digits  for  transmission.   The  error 
will  be  calculated  for  a  decoder  which  makes  an  intermediate 
decision  on  each  bit.   As  noted  this  is  not  the  most  efficient 
use  of  the  channel,  but  it  is  considerably  easier  for  decoding 
and  for  error  calculations. 

4.   Representative  Communications  Systems.   A  simple  system 
which  can  transmit  the  necessary  3°  bits  per  second  is  com- 
posed of  five  parallel  binary  channels.   Each  one  sixth  of  a 
second,  a  five  bit  symbol  is  transmitted,  one  bit  per  channel. 
To  avoid  shipping  noise  a  minimum  frequency  of  900  cps  will 
be  used.   Table  4  lists  the  individual  crossover  probabili- 
ties for  these  five  channels.   The  orthogonal  signal  pair  are 
12  cps  apart  to  allow  a  1/6  second  signal  interval.   The 
channels  are  separated  by  20  cps,  to  permit  easier  tuning, 
in  case  of  a  doppler  shift  which  could  be  as  high  as  15  cps 
for  two  ships  at  15  knots. 

Table  4 
Channel  #  frequencies  ?(c  ) 

1  900  &   912  0.0014 

2  932  &   9^4  0.0019 

3  964  &  9?6  0.0024 

4  996  &  1008  0.0030 

5  1028  &  1040  0.0038 
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The  error  rate,  per  five  bit  symbol  for  this  transmission 
system  can  easily  be  computed  as: 

P(<f)  =  1  -  P(no  bit  errors) 
For  independent  bit  errors  the  error  probability  is  easily 
computed  to  be  0.013  per  five  bit  symbol.   Due  to  receiver 
imperfections,  and  variation  of  the  channel  from  the  model 
which  will  be  discussed  in  the  next  section,  this  can  only 
be  considered  an  order  of  magnitude,   therefore,  methods  of 
improving  the  reliability  will  be  investigated.   Shannon's 
coding  theorem  discussed  in  the  last  section  indicates  that, 
the  error  rate  may  be  improved  without  an  increase  of  power. 

The  easiest  method  to  implement  for  this  channel  is 
not  normally  thought  of  as  coding,  though  it  is  in  reality 
a  simple  code.   This  is  frequency  diversity.   For  illustration 
fifteen  parallel  channels  will  be  used,  with  each  bit  being 
transmitted  over  three  channels  instead  of  one.  '?o   maintain 
the  same  transmitter  power,  each  channel  will  use  only  one 
third  of  the  power  as  for  the  simple  system  binary  channels. 
The  error  probabilities  for  the  fifteen  channels  are  listed 
in  table  5»   The  frequencies  are  spaced  as  for  the  simple 
five  channel  system. 

Table  5 


Channel 

l 

?(<£) 

Channel 

# 

?(£) 

] 

.0025 

9 

.0264 

2 

.0059 

10 

.0324 

3 

.0074 

11 

.0385 

4 

.0090 

12 

.  0446 

5 

.0117 

13 

.0550 

6 

.0149 

14 

.0705 

7 

.018*+ 

15 

.0840 

8 

.0224 

2    / 
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A  block  diagram  of  the  triple  diversity  receiver  is  shown  in 
figure  1^.   The  wA  are  weighting  functions  which  depend  on 
the  signal-noise  ratio  for  the  particular  binary  channel  (33) 


W 


i 


=  Ei/Ni 


For  the  fifteen  channels  being  used,  they  range  from  0.99 
for  channel  number  one  to  0.86  for  channel  number  fifteen. 
For  a  simple  equal  energy  L-fold  diversity  system  the  error 
probability  can  be  computed  from  the  following  formula  (3*0: 

This  is  plotted  for  L  equal  to  1,  2,  and  3  in  figure  15.  As 
an  example,  if  p  is  0.0^,  three-fold  diversity  yields  a  bit 
error  probability  of  0.00065.   For  a  five  bit  symbol  this 
would  become  an  error  rate  of  0.0033*   Assuming  that  p  equal 
to  .0^  is  a  good  average  from  table  5f  this  is  a  four  times 
improvement  over  simple  five  channel  transmission.   A  general 
error  calculation  method  for  unequal  energy  and  noise  on  the 
different  diversity  channels  is  worked  out  in  Appendix  D. 

If  data  is  to  be  transmitted  with  an  accuracy  better 
than  this  diversity  system,  coding  will  be  required.   One 
method,  block  to  block  coding  will  be  discussed  now.   Other 
methods  will  be  compared  to  it  in  the  next  section.   A  prac- 
tical block  to  block  coding  system  is  the  parity  check  code. 
For  illustration,  each  five  bit  symbol  will  be  encoded  into  a 
fifteen  bit  code  word.   The  first  five  bits  of  the  code  word 
are  the  five  information  bits;   the  remainder  are  modulo  two 


31 


channel        / 
receiver 

w, 

m0 

Out  f>u 

t 

w, 

Outp 

A 

channel    *Z 
receiver 

w2 

*y>. 

r  (i )          > 

outpu 

t 

Wz 

n*i 

Outpu 

t 

channel        3 
receiver 

w, 

Output 

w, 

m, 

*-/"; 


output 


Figure  14.   Three-fold  diversity  receiver. 
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Figure  15.   Binary  receiver  error  probabilities  for 
the  Ray lei gh  fading  channel. 
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sums  of  several  of  these  five  bits.   Proper  selection  of  the 
code,  and  use  of  an  appropriate  receiver,  will  permit  cor- 
rection of  any  three  bit  or  lesser  error  in  the  reception 
of  the  fifteen  bit  code  word  (35  )•   F°z*  this  code  the  error 
probability  than  bccoi.-.es  : 

P()  =  l-p(no  errors)  -  p(one  error)  -  p(two  errors)  -  p (three 

errors ) 

Carrying  out  this  calculation  for  the  data  of  table  5t  the 

error  probability  is  approximately  0.004  per  five  bit  symbol, 

or  approximately  the  same  as  for  triple  diversity.   -Further 

improvement  can  be  accomplished  by  coding  by  use  of  a  longer 

block  length,  for  example,  coding  two  five  bit  symbols  into 

a  thirty  bit  code  word,  or  three  symbols  into  a  ^5  kit  word. 

5.   Suitability  of  Coding  'techniques .   In  the  last  section 
simple  transmission,  frequency  diversity  and  a  short  block 
length  parity  check  code  were  compared.   There  was  a  signi- 
ficantly smaller  error  using  diversity  or  the  coding  scheme. 
To  achieve  a  still  smaller  error  probability,  much  larger 
block  length  or  other  types  of  codes  would  be  required. 
Various  coding  techniques  will  now  be  discussed,  which  may 
be  used  if  the  added  complexity  can  be  justified  by  the  im- 
proved performance. 

A  block  to  block  parity  check  code  was  used  as  an 
example  in  the  last  section.   For  this  type  of  coding  system, 
the  number  of  computations  required  by  the  decoder  increases 
roughly  as  the  cube  of  the  block  length  (3&).  and  as  noted, 


34 


only  specific  error  combi nations  can  be  corrected.   This 
limits  its  capability  of  correcting  a  large  burst  of  errors, 
except  for  very  long  codes.   In  general,  block  code  receivers 
make  a  decision  on  each  bit  before  decoding.   This  means  that 
the  confidence  level  in  each  bit  is  lost.   The  performance 
can  be  improved  by  finer  quantization  in  the  receiver  before 
decoding.   It  can  be  shown,  however,  that  three  level  quanti- 
zation will  achieve  half  of  the  possible  improvement,  and 
that  finer  quantization  than  this  gives  an  increasing  small 
improvement  (37 )•   Three  level  quantization  is  best  described 
by  the  binary  erasure  channel  shown  in  figure  16.  For  this 
channel  the  bias  level  must  be  properly  set  so  that  both  p  and 
q  are  small,  where  p  is  the  crossover  probability,  and  q  is 
the  erasure  probability.   This  type  of  channel  is  more  dif- 
ficult to  instrument  than  the  simple  binary  channel  and  has 
not  been  used  as  extensively  (38). 

Some  of  the  drawbacks  of  block  to  block  coding  may 
be  met  with  convolutional  codes  (39).   A  simple  encoder  is 
shown  in  figure  17  (40).   Each  transmission  interval,  v 
parity  checks  on  a  k  bit  register  are  transmitted.   Between 
intervals,  one  bit  is  moved  out  of  the  k  bit  register  on  the 
right,  and  one  new  information  bit  is  moved  in  at  the  left. 
Each  information  bit  is  in  the  transmission  sequence  for  k 
intervals,  and  there  are  an  average  of  v  transmitted  bits  per 
information  bit.   Since  at  least  five  information  bits  must 
be  transmitted  on  fifteen  channels  six  times  per  second,  v  is 
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Figure   16.      The  binary  erasure   channel. 
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Figure  17.  Convolutional  encoder. 
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limited  to  three  for  the  channel  being  described. 

Two  good  decoding  schemes  are  known  for  convolutional 
codes:   sequential  decoding  (41),  and  threshold  decoding  (42). 
In  the  former  an  n  bit  message  is  treated  as  a  tree  with  n 
sequential  nodes.   A  very  good  reliability  can  be  achieved 
if  a  computer  with  a  large  enough  memory  to  backtrack  through 
many  paths  along  the  tree  is  used.   This  may  be  2^  bits  or 
more  (39).   Once  an  error  is  made  because  backtracking  is  no 
longer  possible,  a  large  block  of  errors  may  be  made.   This 
is  satisfactory  in  a  feedback  system  where  repeats  may  be 
requested.   Since  the  transmission  time  over  a  fifty  mile 
hydroacoustic  channel  is  approximately  a  minute,  this  type  of 
channel  is  best  suited  for  broadcast  or  prearranged  time 
sharing,  with  little  or  no  possible  feedback.   A  simpler 
decoding  scheme  for  convolutional  codes  which  avoids  these 
block  errors  at  the  expense  of  a  higher  average  error  rate 
is  threshold  decoding. 

The  threshold  decoder  uses  only  shift  registers,  modulo 
two  adders,  and  a  threshold  device.   For  good  results  each 

bit  must  remain  in  the  register  for  more  than  one  fifteen 
bit  transmission  time,  since  for  this  short  coding  length 
the  results  are  about  the  same  as  for  triple  diversity  as 
noted  on  figure  15  (43).   For  short  block  lengths,  the  results 
with  threshold  decoding  on  a  convolutional  code  are  somewhat 
better  than  for  block  to  block  coding,  with  no  more  complexity. 
This  is  because  it  can  take  advantage  of  the  a  posteriori 
probabilities  in  the  binary  channels. 
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Part  III 
DISCUSSION 

1.  General.   In  this  section  the  assumptions  leading  to 
the  channel  model  and  the  resulting  communications  system 
will  be  discussed  first.   The  importance  of  the  assumptions 
to  the  resulting  communications  system,  and  their  expected 
variations  will  be  brought  out.   Next  the  hydroacoustic  con- 
ditions in  various  ocean  areas  at  different  times  of  the  year 
will  be  discussed. 

2.  Assumptions .   The  most  important  assumptions  made  in  the 
preceding  development  have  been  the  following: 

a.  There  is  no  intersymbol  interference. 

b.  The  raultipath  reception  produces  a  Rayleigh 
amplitude  distribution. 

c.  During  each  signal  interval,  the  received  phase 
and  amplitude  are  constant. 

d.  The  receiver  is  synchronized  to  the  transmission 
interval. 

e.  The  noise  is  white  and  Gaussian. 

f .  Errors  on  different  binary  channels  are  independent. 

Assumptions  a.,  c,  and  d.  are  related.   To  minimize  inter- 
symbol interference,  the  signalling  interval  must  be  long 
enough  for  most  of  the  energy  transmitted  during  a  given 
interval  to  be  received  during  an  equal  delayed  interval  at 
the  receiver.   The  time  used  is  based  on  the  calculation  in 
appendix  G.   The  intersymbol  interference  can  be  treated  as 
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an  additional  noise  component  which  increases  the  crossover 
probability,  or  if  it  can  be  estimated,  used  in  a  receiver 
with  memory.   By  making  the  signalling  interval  long  enough 
for  intersymbol  interference  to  be  negligible,  the  interval 
may  become  too  long  to  assume  a  constant  phase  and  amplitude 
during  the  interval.   This  problem  can  be  corrected  by  using 
time  diversity  reception  during  each  signal  interval.   Since 
a  small  amount  of  intersymbol  interference  is  unavoidable, 
and  also  because  the  transmission  time  between  the  terminals 
will  vary,  it  will  be  impossible  to  obtain  perfect  synchroni- 
zation at  the  receiver.   Range  changes  may  be  corrected  by 
a  doppler  circuit.   Transmission  path  variations  should  be 
fairly  small  so  that  the  time  variations  due  to  this  phenomena 
wi  1  be  small  compared  to  the  signal  interval.   In  this  case 
these  variations  are  also  correctable.   In  summary,  these 
assumptions  can  be  met  as  well  as  desired  by  proper  selection 
of  the  signalling  interval  and  the  use  of  time  diversity 
reception  if  required. 

Experimental  data  indicates  that  the  Rayleigh  ampli- 
tude distribution  is  a  good  assumption  much  of  the  time,  but 
that  it  may  differ  (19f20).   For  long  ranges  these  differences 
are  most  evident  in  the  higher  moments,  and  do  not  differ 
significantly  in  the  mean  or  mean  square.   At  short  ranges 
only  a  few  paths  may  be  of  interest,  and  in  this  case  the 
distribution  is  more  favorable  than  a  Rayleigh  distribution 
(19).   Figure  18  compares  a  typical  short  range  amplitude 
distribution  with  the  Rayleigh  distribution. 
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Figure  18.   Typical  amplitude  distributions. 


a.   Upslope  from  the  transmitter 


b.  Downslope  from  the  transmitter. 


Figure  19.   Effect  of  a  sloping  bottom  on  communications 
via  the  bottom  reflection  channel  in  shallow  water. 
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Although  most  ambient  sea  noise  can  be  considered 
white  and  Gaussian  for  the  narrow  band  binary  channels,  other 
forms  of  noise  may  be  present.   Rotating  machinery  harmonics, 
6r  marine  life  may  affect  only  one  or  two  frequencies.   Kost 
of  these  noises  are  most  severe  under  1000  cps,  but  a  source 
on  a  receiving  ship  is  potentially  dangerous.   A  careful  noise 
analysis  of  potential  receiving  vehicles  will  permit  optimum 
frequency  selection  with  respect  to  the  noises.   Impulse 
noises  may  also  be  present.   These  can  be  minimized  by  clip- 
ping the  signal  at  a  level  which  depends  on  the  average  im- 
pulse noise  power  and  the  expected  signal  levels. 

In  arriving  at  L   system  error  probabilities,  errors 
were  considered  independent  of  other  channels.   If  these 
binary  errors  are  not  independent,  then  the  resulting  system 
error  probabilities  x\Till  be  different.   If  for  example  the 
errors  on  two  of  three  diversity  receivers  are  correlated, 
the  resulting  bit  error  is  only  about  as  good  as  double  and 
not  triple  diversity.  No  good  data  is  available  on  how  closely 
correlated  transmission  over  adjacent  frequencies  is. 

3.   Ocean  Conditions  and  the  Channel  Model.   Several  factors 
affect  the  prevailing  ocean  hydroacoustic  conditions  (40). 
The  temperature,  and  hence  the  vertical  profile  of  the  velo- 
city of  sound  vary  with  both  latitude  and  time  of  year,  and 
in  some  cases  the  time  of  the  day.   The  depth  of  the  water 
and  the  bottom  profile  affect  bottom  reflection  transmission. 
Weather  conditions  affect  both  noise  and  surface  reflection 
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losses.   All  the  calculations  made  have  been  for  a  moderate, 
state  three,  sea.   Heavier  weather  conditions  reduce  the 
reliability  considerably. 

Surface  trapped  channels  occur  only  when  the  water 
temperature  is  isothermal  or  perhaps  increasing  with  depth, 
'■'•'his  condition  is  somewhat  unstable  and  can  occur  only  when 
the  surface  waters  are  being  cooled.   Experimental  data 
indicates  that  these  conditions  exist  most  of  the  year  in 
polar  regions,  during  winter  in  mid-latitudes,  and  only 
rarely  in  tropical  waters.   In  mid-latitudes  the  layer  depth 
will  be  from  1000  feet  to  1500  feet  in  the  winter,  but  a 
sharp  thermocline  will  exist  near  the  surface  in  the  summer. 
Since  water  depth  is  not  important  to  the  surface  channel, 
the  model  developed  will  apply  to  all  high  latitude  ocean 
areas  in  the  winter,  but  only  to  shallow  water,  or  polar 
waters  during  the  summer  months. 

The  bottom  reflection  channel  can  also  be  quite 
variable  as  to  results.   For  any  given  depth,  the  bottom 
slope  is  quite  important,  as  is  illustrated  in  figure  19. 
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Part  IV 

CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Conclusions .   In  Part  I  it  was  shown  that  the  signal- 
noise  ratio  for  the  surface  channel  and  the  shallow  water 
bottom  reflection  channel  could  be  estimated.   It  was  further 
shown  that  the  appropriate  model  was  a  Gaussian  random  pro- 
cess, in  which  the  received  signal  had  a  Rayleigh  amplitude 
distribution  and  unknown  phase  with  additive  white  Gaussian 
noise.  Using  this  model  an  appropriate  binary  receiver  was 
found  to  be  a  quadrature  correlation  device.   A  combination 
of  five  of  these  receivers  could  be  used  to  send  data  at  60 
wpm  with  a  moderate  error  rate.   Improvement  in  the  error 
rate  was  best  obtained  by  frequency  diversity,  and  could  be 
further  improved  by  coding  at  the  expense  of  complexity  and 
cost.   In  Part  III,  it  was  shown  that  the  model  was  reasonable 
for  shallow  water  any  time  if  the  bottom  did  not  slope  away 
from  the  transmitter  too  fast,  for  mid-latitudes  in  the 
winter,  and  for  polar  regions  most  of  the  year,   '^hese 
results  are  applicable  for  moderate  weather  conditions,  and 
ranges  of  up  to  fifty  miles. 

2.  Recommendations .   Although  the  basic  form  of  the  receiver 
could  be  developed,  the  actual  parameters  could  only  be 
estimated.   A  better  understanding  of  these  parameters  is 
needed  for  the  actual  development  of  a  suitable  system. 

Some  items  of  interest  which  warrant  additional  experimental 
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research  are  the  following: 

a.  the  noise  spectrum  of  potential  receiving 
vehicles . 

b.  the  range  and  distribution  of  multipath  time 
delays  for  various  channels. 

c.  the  length  of  time  during  which  phase  and  ampli- 
tude of  a  multipath  signal  will  remain  coherent. 

d.  the  correlation  of  errors  on  binary  channels  of 
adjacent  frequency. 
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APPENDIX  A 

1.   Surface  Channel  Sample  Calculations. 

For  these  calculations  the  following  assumptions  will  be  made 

Isothermal  water  with  a  velocity  gradient  of  0.C052  m./ 

s  e  c . /m . 

Frequence  ...  1  kc . 

Range  . . .  50  miles. 

Layer  depth  ...  30.5  meters. 

Sea  state  ...  3 

Source  ...  185  db  (ref.  0.00002  newt./m, )  at  lA  m. 

radius.   This  is  173  db  at  a  1  meter  radius. 

a.  H  =  — c, —   (sin  G  _  gin  Qj 

a  cos  0. 

0=0 

c0  =  1500  m/sec  . 

a  =  0.0052  m/sec/m. 

Y  =  30.5  m 

z 

cos  0o  =   c  ^      1-Q. 

c.  +aY  ~  2 

0o  =  0.0152  radians  =0.37  degrees 
R  =^380  meters 

b.  Spherical  spreading  loss. 

L  =20  lore  3/2  =  2  0  log,  2190  =  67  db. 

c.  Cylindrical  spreading  loss. 

L  =  10  log  50  miles  =  10  log  36.7  =  15-7  db. 
"   R/2 

d.  Surface  reflection  loss. 

Number  reflections  =  50  miles  =  9.2 

2R 


i±5 

Number  reflections  =10 

Loss  per  bounce  =  0.95  db  for  a  4  foot  wave  height 

L  =  0.95  X  10  =  9.5  db. 

e.  absorption.   Negligible  for  this  range  and  frequency. 

f.  Total  losses.   6?  +  15.?  +  9.5  =  92  do. 

g.  Signal  level.   1?3  -  92  =  81  db  (ref.  0.00002  newt./m.) 

2 .   Bottom  Reflection  Channel  Calculations. 

The  same  assumptions  are  made  as  for  the  surface  channel 

except  there  is  no  layer,  and: 

Water  depth  ...  300  meters,  flat  bottom. 

Velocity  gradient  ...  0.0?5  m./sec/m. 


a.   R  = 


(sin  0o) 


R  =  3500  meters 

b.  Spherical  loss. 

L  =   20    log/fl3500   =   71   db. 

c.  Cylindrical  loss. 

L  =  10  log  50  miles  =  10  log;<)23  =  13.6  db. 
"      3500 

d.  Bottom  Loss. 

Number  reflections  =  50  miles  =  11.5 

2R 

Number  reflections  =  12 

L=  12  X  0.3  =  3-6  db. 

e.  Surface  loss. 

L  =  12  X  0.95  =  11.4  db. 

f.  Total  loss.   71  +  13.6  4-  3.6  +  11.4  =  100  db. 

g.  Signal  level  .  173-100  =  73  db.  (ref.  0.00002  newt./m.) 


k6 

APPENDIX      B 

1 .  Cavitation, 

Cavitation  occurs  when  the  acoustic  pressure  plus  the  ambient 
pressure  is  less  than  the  vapor  pressure  of  water.   The 
following;  numerical  values  will  be  used  in  this  calculation: 

vapor  pressure  =  0.6  psi  at  85  degrees  P.  (decreases  with  T) 

atmospheric  pressure  ...  1^.7  psi 

water  depth  =  20  feet 

a.  Static  pressure. 

p  =  1^.7  +  0.W*  X  2  0  =  23.5  psi 

b.  Allowable  acoustic  pressure. 

p  =   23.5    -   0.6   =  22.9   psi  =   1.58  X  10  newt./m. 

c.  Sound  pressure  level. 
SPL  =20  log„p/p(ref ) 
p(ref)   =   0.00002   newt./nu 

SPL  =   20   log,/?. 9  X   10*  =   198  db. 

Allowing  3  db  for  peak  versus  mean  square,  the  SPL  is 

limited  to  195  <lb  (ref.  0.00002  newt./m*. ) 

2 .  Non- Linearity. 

For  the  acoustic  approximation  we  require: 

u<)u  ~  hi 
i  r   <)  t 

Assuming  spherical  harmonic  waves,  the  following  equations 

apply: 

P  =  2»    e;)  (".t-kr) 
r 


U  =pH-ik)  Pd 

tr  J*7  -yre. 


j(*vt-kr) 


w 


Taking  the  indicated  derivatives  we  now  obtain: 


p 


1  +  jk  +  k      «    j 


*. 


jft  [rv     r 


let    f ,  =  -gfr    =    1000   cps 


r   =   1  meter 


u). 


then,  since  k  -  -~  >  the  bracketed  term  on  the  left  is  less 
than  10.   1he  inequality  can  be  simplified  to 

2-    «  Jt 

e. 

In  sea  water,    £,  «  1000  kg/m. 

o>.  =  1*  77  X  10   at  f#  =  1  kc. 

so  that  the  inequality  becomes: 

P  «  k  w*  X  10 
8 

Letting  p  =  10  ,  we  obtain  the  following  sound  pressure  level 
SPL  =  20  log,  p/p(ref ) 


=20  log  10 

10 


0.00002 


=  25^  db  (ref.  0.00002  newt./m.) 
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APPENDIX  C 

1.  xi~re   Dispersion.   The  maximum  difference  in  the  length 
of  time  it  takes  a  signal  to  reach  the  receiver  via  two  dif- 
ferent multipaths  limits  the  rate  at  which  signals  may  be  sent. 
This  difference  is  due  to  sound  velocity  fluctuations,  and 
path  length  fluctuations.   This  time  difference  will  now  he 
estimated  for  the  special  case  of  a  bottom  reflecting  channel. 
Variations  will  be  less  severe  for  surface  channels  or  for 
shallower  water.   It  will  be  assumed  that  the  average  velocity 
of  sound  differs  by  a  maximum  of  four  meters  per  second  over 
different  paths.   This  corresponds  to  an  average  temperature 
difference  of  1.0°F.   This  is  a  reasonable  value  since  va- 
riations over  several  miles  are  normally  quite  small  in  open 
waters.   It  will  further  be  assumed  that  the  ray  paths  are 
straight  lines.   This  is  illustrated  in  figure  20.   Since  0 
is  limited  to  10°  -  15°  at  the  bottom  for  good  reflection 
coefficients,  an  average  value  of  6°  will  be  chosen.   This 
is  reasonable  since  the  most  significant  paths  are  those 
leaving  the  transducer  within  a  few  degrees  of  the  horizontal. 
All  these  paths  will  reach  the  bottom  at  approximately  the 
same  angle,  so  that  the  average  variation  will  be  only  about  1°. 

Using  these  values,  a  range  of  50  miles,  and  a  sound 
velocity  of  1500  m/sec,  the  time  variation  will  now  be 
computed. 

T  =  D/V  =  slant  range/velocity 

R  =  D  cos  0  =  horizontal  range 
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Figure   20.     Model   for  time  dispersion  calculation 
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dT  =    -  D  dv   _  R      tan   Qd0 
vA  v 

^  -  1  dv   -  R   Qd© 
v*  v 

inhere 

R  =  50  miles  =  80,500  meters 

V  =  1500  m/sec. 

dV  =  *Mn/sec. 

0  =  ,104  radians 

dO  =  .0174  radians 

Carrying  out  the  calculation  the  time  variation  is  0.2^0  sec. 
This  is  a  variation  from  the  straight  line  path,  however. 
Even  the  shortest  paths  will  range  from  the  straight  line  by 
,0k   -    .06  sec.   This  leaves  a  total  variation  of  about  0.2 
seconds  which  compares  favorably  with  the  only  experimental 
data  available  (23). 
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APPENDIX  D 

1.   Diversity  5rror  Calculation.   In  part  II  it  is  noted 
that  the  diversity  receiver  makes  a  decision  that  me  is 
sent  if: 

2L  2L 


where  the  a..l      are  statistically  independent  zero  mean  Gaussian 
random  variables.   Since  the  VJji    are  all  approximately  equal 
they  will  be  neglected  and  the  error  probability  will  be 
computed  for  making  the  decision  on: 

PL  2i 

?  «5      >     i  - 

j  =  /  £*i 

Ey  pairwise  addition  of  the  quadrature  components  on  each 

channel,  the  decision  can  also  be  represented  as: 

l$    >  it 

where  the  b-  and  the  b^ are  statistically  independent  Rayleigh 
distributed  random  variables,  with  the  following  probability 
distributions . 


P(bj  )  =  fo  exp  _  2b^_  o   <  bj  < 


oO 


p(b,-  )   =  Lj±i_  eXD   _  2^  0  £   hi  <    °° 

The   error   probability  can  be   easily  calculated   if  the   pro- 
babilities   of   the   two  sums      are   known.      These  will  noi\T   be 
put   in  a  closed   form.      Define  a  new  variable: 

then   it    is    easily  shown   that'. 

P<*j  )    =        2        ex.Q   _  2*i__ 
Ej  +Nj  *         Ej  +N,- 


- 
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The  probability  distribution  of  the  sum  of  the  k.-   can  be 
found  by  noting  that  the  characteristic  function  of  the  sum 
of  independent  random  variables  is  the  product  of  their  in- 
dividual characteristic  functions.   The  characteristic  fun- 
ction of  k  is  : 

fp(kj)  e-skJdkj=:  2 (  1     . 

Therefore  the  characteristic  function  of  the  sum  of  the  k.-   is 

fptfk,-  R,)  e-sR,  dH  =  jj-(     2    }  (    1 

J*  /=/  Ej  +NJ      s+2/E4  +N,-  ; 

By  a  partial  fraction  expansion  it  is  easily  shown  that  the 
probability  density  function  of  the  sum  of  the  k;   is  the 
sum  of  L  exponentials. 

Knowing  the  probability  functions,  the  probability  of 
error  can  easily  be  calculated.   If  ui0  is  sent,  and 

d=l 


L 


*jt  =  R/ 


then  the   error   probability   is: 
p(£)   =      (p(R,  )     fpttO    dR0dR, 

This  integral  can  be  solved  in  general.   For  the  special  case 
of  two-fold  diversity  it  takes  the  following  form: 

p(<£ )  =  1  -  /  K,Kv\f/  K,  \  /__!_   -  .  _1\  +  /K, V 1_  - 


Kj-Ky/AKy  -K2/lKa+KY      Kz+K,/         \  K,  -Kjl  K,  -rZ£3      K,  +K¥// 
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where : 

K; 

= 

2 

S,  +N, 

-z 

= 

2 

E2  +N2 

A3 

= 

2 
N, 

K 

= 

2 
N2 

j-'/.is  last  result  is  restricted  to  cases  where  K,  ^  K2and 
Kj^K^   If  one  or  both  of  these  pairs  were  equal,  the  result 
would  be  found  in  a  similar  way,  but  would  be  slightly  dif- 
ferent.  Figure  21  shows  the  effect  of  varying  signal-noise 
ratio  on  two -fold  diversity. 


s, 


■k 


looo 


3/5"-- 


100  -- 


■31.5- 
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3.2-- 


K  = 


.  Vm 


V/V, 
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Figure  21.  Diversity  error  with  unequal  signal-noise  ratios 
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